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There is an increasing tendency towards the use of a lesser number of controllers in 
industry. This is so as the inclusion of an economical optimization layer becomes 
simpler. However a reduced number of controllers means that they are large ones. 
These must be adequately tuned, which can be accomplished through their state space 
representation. Several ideas for obtaining a state space representation of the class of model 
predictive controllers (MPC) have appeared in the literature like the approach of Li et al\ 
Nevertheless, most of them lead to a system with thousands of states if the controller becomes 
large. As a consequence, the computation of eigenvalues or singular values to evaluate the 
stability of the system becomes prohibitive. In this paper a new state space representation of the 
MPC controllers is introduced based on the continuous step response of the system given by a 
transfer function. The idea is to represent the step response in a parametric form. The number 
of states of this representation is of the same order as the identified transfer function however 
large the optimization horizon of the chosen predictive controller might be. The method 
includes integrating and open loop unstable processes. Besides that the model time delay 
appears explicitly and is independent of other parameters. The computation work involved in 
the closed loop analysis is thus minimized. Also, this new representation can be used to 
robustly tune the MPC controllers without the use of extensive simulations. 

Keywords: model predictive controllers; robust control 



I. INTRODUCTION 

For almost two decades, linear model predictive controllers 
have been extensively used in industry (Richatet et ai 2 t 
Cutler and Ramaker 3 , Prett and Gillette 4 ). The main reason for 
their acceptance is their capability of enforcing constraints on 
variables and their ability to capture the economic benefits by 
an optimal operation of the plant. However, two additional 
interesting properties of this type of controller should not be 
forgotten, namely, their stabilizing properties and their 
multivariable character. 

Both of these are important for achieving optimal 
operation of a plant. In recent times the demand for plant 
optimization has grown enormously. The plant economic 
objectives are more easily defined for the full unit. At the 
same time it is easier to implement the optimization 
commands (normally given as set-points to the controllers) 
if there are a small number of controllers. The ideal situation 
would be to have a single controller for the whole plant. If 
one considers a real process, e.g. if the optimization of a 
crude distillation unit or a fluid catalytic cracking unit is to 
be undertaken, the number of manipulated inputs may 
amount to fifty and the number of controlled outputs may be 
near a hundred. From the implementation point of view 
there are also advantages in increasing the size of the 
controllers. The controller commissioning gains more 
flexibility since this work can be done by switching on the 
controller per section of the plant. 

Thus it can be seen that being able to deal with a large 
controller is what should be aimed for. On the other hand, 



most of the industrial MPC packages such as DMC (Cutler 
and Ramaker 3 ), MPHC (Richalet et a/. 2 ), MAC (Rouhani 
and Mehra 5 ) and STAR (Guettler et aL 6 ) still use the same 
concepts of the early predictive controllers, i.e. they are 
based on the step response model of the process. Therefore, 
the number of coefficients in the models used by the 
controller may be greater than 100. This is necessary if a 
good performance of the system is required. 

The MPC packages when applied to a large number of 
inputs and outputs may have to deal with some thousands of 
parameters and model coefficients. With the available 
computer power and optimization software, this might not 
be a major problem for the calculation of the control actions. 
However, as far as the tuning of the controller is concerned, 
problems may arise. The work involved in the usual practice 
of simulating the system for typical disturbances and model 
uncertainties tends to become overwhelming, as the number 
of possible situations will increase at an exponential rate. 
Frequently, analysis tools cannot be found for the design 
procedure of the MPC. The availability of these can, 
however, simplify this procedure. 

Li et al} showed how to represent the usual uncon- 
strained MPC strategy in the state space representation. By 
doing so the way to the closed loop analysis of the predictive 
controllers was open. This analysis encompasses many 
features, as for example: 

• how to evaluate the robustness of the controller with- 
out extensive simulation it is possible to evaluate the 
performance of the controller submitted to model mismatch 
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or how to select its parameters so that the performance 
remains acceptable for a wide family of identified models. 
• how to evaluate whether the controller maintains a 
suitable response, e.g. by constraining the closed loop poles 
inside a closed region of the complex plane. 

So, if the system is in the state space form, an extensive 
library of methods and results related to the linear system 
theory could be incorporated as analysis tools for the model 
predictive (MP) controllers. The available procedure (Li 
et al} and Lundstrom et al?), nonetheless, produces for a 
large system a model with thousands of states as each step 
coefficient corresponds to a new state. Under this scenario 
there is a clear incentive to obtain a more economic 
representation of the system! 

In this paper a more compact state space representation of 
the class of MP controllers (named ROSSMPC) that uses 
step response based internal models is presented. In section 
II the existing approaches is briefly reviewed and it is shown 
why they may become of little use for an extensive closed 
loop robust analysis. In section III a state space representa- 
tion is presented and the way the order of the model can be 
reduced is shown. In section IV some examples are 
presented where it is shown how ROSSMPC can be used 
in the analysis and design of MPC controllers. For example 
it will be shown that the DMC controller can have the 
same performance as the OBMPC controller if properly 
designed, contrary to what was claimed in Lundstrom et al. . 
In section IV it is also shown that a subsystem of the well 
known Shell problem (Prett and Garcia l0 ) can be stabilized 
with the present ROSSMPC formulation with a much lesser 
optimization horizon than the ones obtained with the 
conventional representations presented in the literature. 

As the controller must retain its good performance in the 
real environment, it is further shown how a new analysis 
tool (Hudchison and Pritchard 11 and Li et al. 12 ) can be used 
to quantify the robustness of the chosen controller. This tool 
is generalized for some structured systems. This is shown in 
section V where this method is introduced as an additional 
tool that can be used together with the existing techniques 
for the design of MPC. In section VI some examples are 
presented of how these tools can be used and the paper 
finally concludes in section VII. 

II. THE EXISTING STATE SPACE 
REPRESENTATIONS AND THEIR LIMITATIONS 

Most of the original MP controllers make use of step or 
pulse response models to predict the system output. The 
coefficients of these models can be obtained by their fitting 
to plant data (Bailey 13 ). The model so obtained is usually of 
a very high order. Time delay is accommodated by placing 
zeros in the first model coefficients. 

One point that must be emphasized is how the time delay 
is enabled to vary in this model. This is due to the fact that in 
the closed loop analysis the effect of the time delay on the 
stability or performance of the system must be considered. If 
it appears interrelated with other model parameters, the 
analysis becomes more complicated. In the step response 
representation changes, in the time delay can be considered 
in the following two ways: 

(a) The new time delay is equal to the old time delay plus or 
minus a multiple of the sampling time. In this case the new 

4 



model is obtained simply by adding or removing leading 
zeros from the old mode). The remaining model coefficients 
are not modified. Therefore the time delay is an explicit 
parameter. 

(b) The new time delay affects the whole response. In this 
case not only the leading zeros are affected but also the 
coefficients of the step or pulse response model. Conse- 
quently it can be considered that in this model formulation, 
the time delay is not an independent or explicit parameter of 
the model. 

Not only is the time delay an important parameter that 
should be made explicit whenever possible, there are other 
intuitive parameters of the system that should be made 
explicit. For example the control engineer has a good idea of 
how much the gains and time constants of the real process 
can change in the real system. This notion is very important 
in the closed loop analysis. Though it is not simple to 
associate these instincts with the coefficients in the step or 
pulse response models. 

With the introduction of the Matlab Identification 
Toolbox® another model became popular, namely the 
pulse transfer function (Ljung 14 ). This is usually a low 
order parametric model where all the intuitive aforemen- 
tioned parameters always appear explicitly except the time 
delay. 

Lee and Yu 19 and Oliveira and Biegler 20 propose a 
general state space structure to represent and analyse the MP 
controllers. In their representation, time delays are not 
considered and they give no indications on how to deal with 
the situation when this parameter can vary continuously. 
Another limitation of their approach is that the resulting 
relation between the predictions of the system output and 
the model parameters is nonlinear. This last aspect 
complicates the robustness analysis of the closed loop. 

Most chemical engineering processes are continuous and 
can be adequately represented by continuous transfer 
functions. For MIMO systems the relation between any 
controlled output v ( £ //?"> and any manipulated input 
u j G IR n " may be represented by equation (1). 



where 



where 






is the process continuous transfer function 


*J 


are indexes that are related to the controlled 




and manipulated variables 


s 


is the Laplace domain 




is the time delay appearing in an explicit 




manner 


Bij and A u 


are polynomials given by: 


B u (s) = 


+ b iJA s + b^s 2 ... + b^ nbt / b - 


Avis) = 


1 + a iJA s + a iJm7 r... + <ii J .na iJ s naiJ 


bij.k and 


^i.j.k are the coefficients of the model 



nbf j and na^ are the orders of the polynomials 



The polynomial Aq(s) is usually assumed to have the 
single roots: r^j, r,j 2 . r iJjul ... This is not a serious 
restrictive assumption, since poles of high multiplicity are 
not frequently found in the process industry. The roots of the 
polynomial A,,/j) correspond to the inverses of the time 
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constants of the system. The process gain and the 
parameters of eventual inverse responses are directly 
associated with the coefficients of the numerator of the 
transfer function, namely 

So it can be seen that the model of equation (1) seems to 
be an adequate model that could be used in a closed loop 
analysis of MP controllers. As the implementation is 
normally carried with the convolution model, the relation- 
ship between these two formulations must be obtained, and 
this is presented next. 

If a step is applied to the system described by equation 
(I), the corresponding step response model can be obtained 
for any instant as given in equation (2). 

and the sampling instant is such that nT > The 
coefficients c iJg are the residues of G^/s and/-^ is given by: 

for g = 0 : f w {i) = 1 



for g * 0 : 



f *W-*/> : r u . * 0 

Q n (3) 

From the above equations it is clear that the step response 
coefficients are explicit functions of the parameters of the 
transfer function model. The gain of the system relating 
input iij and output^ is related linearly with the coefficients 
c iJtg . So if there is any uncertainty in the parameters of the 
model (due to e.g. plant mismatch) its influence can be 
easily evaluated. For instance, an increase of a given 
percentage in the plant gain corresponds to the same 
percentage increase in c iJg for all g. 



The Closed Loop Formulation of OBMPC 

For open loop stable processes, the state space model 
presented by Li et al.\ allows the state predictor of MPC to 
be written as: 



where 

Au k = iu - if.. 



(4) 
(5) 

(6) 



M ■ 



0 


7 -r 


0 . 


0 


0 " 




" s z - 


0 


0 




0 


0 






0 


0 


o ■ 


• **, 


0 


s = 




0 


0 


0 


0 


J„ 






0 


0 


0 


0 









(7) 

S, is the i th coefficient of the step response 

n h is the number of coefficients of the step reponse model 

R is the plant output correction gain. It can be easily 



calculated following the procedure of the Kalman 
filter . In this case Lundstrom et al. 1 designate the 
controller OBMPC (Observer Based MPC) 
y k is the real measured output vector at instant k 
N is a matrix which collects the first n y components of the 
state vector 



N = 



A, 0 ... 0 



If a model with the same structure of the predictor but with 
different parameters is used to represent the real plant in the 
closed loop analysis, a total of 2n v .n h stales will be 
necessary. Since n h is usually about 100 and « v may go up 
to 50 for a large controller, the method would generate near 
ten thousand states. This is certainly not suitable for 
practical purposes. 



III. THE REDUCED ORDER STATE SPACE MODEL 
PREDICTIVE CONTROLLER (ROSSMPC) 

In the previous section it was shown how a state space 
representation of the predictive controllers based on step 
responses curves can be obtained. The approach shown 
there was that of Li et al. 1 It could be seen that the number of 
states was very large. Moreover, it was equal to the order of 
the convolution model, i.e. equal to the number of 
coefficients. So the number of states does not depend on 
two interesting tuning parameters of the controller, namely 
the control and optimization horizons. One attempt to 
reduce the number of states was made by Hovd et al. 9 
and Lundstrom et al. 7 The authors have noticed that a 
truncation of the step response curve will reduce the number 
of states. In order to preserve the information of the 
truncated part, a parametric first order model was chosen to 
approximate the tail of the curve. This approach, besides 
being approximate and consequently not useful for closed 
loop analysis, is only valid for open loop stable processes. 

In regard to the aspects considered in the last paragraph, 
the points that must be taken into account for obtaining a 
good representation of the controller can be summarized as 
follows: 

• as the optimization horizon can be much smaller than the 
order of the convolution or step response model, it would be 
desirable to have a representation that would be dependent 
on the optimization horizon. 

• the formulation of the controller should be valid for both 
stable and unstable open-loop processes. 

In addition to these aspects, another point can be outlined. 
Rowlings and Muske 16 showed that under the hypothesis of 
ideal model and open loop stable systems, the MPC is 
intrinsically stable for any tuning parameters if the 
optimization horizon tends to infinity. So it seems reason- 
able that for the real process where there is model 
uncertainty, a long prediction horizon tends to stabilize 
the system. This justifies a practical rule that recommends 
the use of an optimization horizon equal to the open loop 
stabilizing time of the system. With this rule there are 
processes for which the MPC requires a large control 
horizon. That means that even if the number of states were 
reduced according to the criterion based on the optimization 
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horizon, this would be not acceptable. So another desirable 
characteristic of the representation would be: 
• the capability of reducing the number of states when the 
optimization horizon must be large for stabilizing purposes 
Our approach, named ROSSMPC (reduced order state 
space model predictive controller) incorporates these three 
desired characteristics. In the rest of this section the 
ROSSMPC formulation will be introduced. This will start 
with the prediction model that is dependent on the 
optimization horizon rather than on the convolution order. 
So it will be shown how to transform the dependency on the 
convolution order into a dependency on the optimization 
horizon. As this will be done it will also be shown that the 
approach is also valid for unstable processes. After this the 
closed loop response is derived together with the control 
action. Finally it will be shown how to reduce the obtained 
representation making use of different fictitious sampling 
periods. 

Building the Prediction Model of the ROSSMPC 
The starting point will be a function whose output is the 
prediction curve of the controlled variables given by a 
weighted composition of the step responses in the inputs. It 
will be obtained on all available information up to the 
instant k. So the ROSSMPC prediction model has the 
general form given by: 

n„ rut,J 

where, /=1 n y J=\ £=1 na tJ , p t and p iJtg are 

the new states of the model that will be called artificial states 
in the remaining of the paper and n is taken such that 
nT > max, ^. In multivariate form equation (8) becomes 



where, 



(9) 



P = 



where 



(12) 



,f r - =0 

Or alternatively in compact multivariate form as: 

1,* = IPhk + OA'* 

IPUuk = Wit* + CD ^ 
where. 



(13) 
(14) 
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D = diag(d UA d u2 ••■ d \A^ [A d ixi d \xi d U2^ l2 

•■• <A.n..1 ^1.".^,^ *" ^2.1.1 ••- ^Ji.-Mrw 

--- d n y A.l d n y A.2 ••• d n y ^t u jui„ yJtM ) 

As low order transfer function models are normally used 
to represent the system, the number of states generated by 
equations (13) and (14) are small. This means that the 
reduction in the number of states when considering 
equations (13) and (14) instead of the states corresponding 
to the coefficients of the step response ranging from the 
optimization horizon (n,) up to the stabilizing horizon (n h ) is 
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In equation (10) ^ is calculated in the instant (nT). 

Next, how the artificial states of this model can be 
obtained will be discussed. Note that equation (8) is very 
similar in form equation (2). This enables the artificial states 
to be updated as follows: 



(10) 



[/>,1*+1/A = [Alt 



(ID 



.1.2 fn y Ajta. y .i /i„.2.l "" A- 2 -'""- *"* A,ji fc .l **' fa,******. 

substantial. An additional advantage of this model repre- 
sentation is that the prediction of the system output can be 
arbitrarily extrapolated including time instants far into the 
future without significantly increasing the order of the state 
space model that represents the closed loop. 

It is not a common feature of the MPC packages found in 
industry to consider not equally spaced prediction time 
instants. However, with the proposed model structure a 
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controller can be designed with a large optimization time 
horizon and few prediction instants which results in a small 
closed loop state space model. This property will be 
commented on later. 

Equation (4) can also be used as the prediction model of 
ROSSMPC. The only difference lies in the definition of the 
matrices M, S and the vector of states. This last is redefined 
so as to incorporate the artificial states of the model in P and 
P. The new matrices M and S are obtained by combining 
equations (4), (13) and(14). Thus: 



OWi/* = Mly]m + SAu k 



(4) 



with 



0 l Hf 0 ... 0 0 
0 0 /„ ... 0 0 



0 0 0 
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(15) 



Equation (4) is used to predict what will happen if no 
other control actions are taken. As the process is to be 
stabilized, it must be known what the states will be like 
when steady-state is reached. Here lies the main difference 
between our approach and others. For the conventional 
MPC strategies (e.g. DMC, MAC, etc) there is the need to 
predict n h future instants, where n h is the time needed for the 
process to stabilize. For real processes this number is 
usually very large (1 00 or larger). In the ROSSMPC strategy 
presented so far there is only the need to predict n,n y + n p 
future values, where n r is the optimization horizon and n p is 
the number of artificial states in P and P and can be obtained 
as: 



"/> = "> + £X>4v 



=1 j=\ 

For real processes the required number n r may be as low 
as 10 and a typical value of na u is 2. Also, for a real process 
n h is of the order of 100. Using the ROSSMPC representa- 
tion for a SISO system, only need 13 states would therefore 
be needed in equation (4), while the conventional state- 
space representations of the MPCs will need n h states, i.e. 
about 100 states. So it can be seen that there is a significant 
reduction in the number of states required. 

Before going further, comments are needed on some 
characteristics of our model. Observing equation (15) and 
the definition of matrix ^ (equation (10)) it can be seen that 
all the terms in matrices M and S preserve a linear relation 
with each of the parameters of the original step response 



model (equation (2)). This is important if the uncertainties 
in the model are to be taken into account. Also their 
continuous form is preserved and this makes it possible to 
search for the worst case condition where the controller 
might have performance or stability problems as will be 
seen later in this paper. 

Again, it is desirable to correct the prediction with 
measured data which is accomplished by the use of equation 
(5). Combining equations (4) and (5) the prediction model 
given in equation (16) is finally achieved. 

Lv] 4+ ,/* +l = (/ - RN)M[y] k/k + Ry k + </ - RN)SAu k 

(16) 

The control actions are calculated so as to minimize an 
objective function. MP controllers assume that equation 

(17) is the objective function. In equation (17) the output is 
predicted in n r future instants and m control actions are used 
to drive the process to the desired states. The effect of the 
future control actions are predicted according to equation 

(18) . 

J = ?i n {H rCvV ~ kXJll 2 + IIAAt/J 2 } (17) 



Mm = IvL* + SAU k 
where, 



(18) 



AJT are diagonal weighting matrices 

yp ^ /K«,™ y +n P tne desired values for the states and 

will be commented on later. 
AU k are the future calculated control actions 

given by: 

i r 



0 



0 
0 

0 



and Si are the same as in equation (15). 

The well known least-squares solution of equation (17) 
with the prediction (18) is presented in equation (19). 

AU k = {S T r T rS+A T A}- l S T r T r{y" - [yU} (19) 

As only the first control action is of interest, the following 
can be we obtained: 

«W= K 0 ... Q]AU k 

So equation (19) can be written in a compact form as 
equation (20). 

Au k = K MPC {y sp - [y] w ) (20) 

A natural question that can arise is how to evaluate y sp for 
the artificial states in P and P. To answer this question, note 
that it is necessary to observe that if equation (4) is to be 
valid at steady state, P has to be driven to zero and P to 
the system output. Consequently, the components of y sp 
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corresponding to P in equation (20) are set equal to zero 
while the components of y sp corresponding to P , have the 
same values as the set points of the outputs. 

The Closed Loop Response of the ROSSMPC 

Taking advantage from the fact that the last n p ~n y terms 
in y sp must be zero and the particular structure of matrix Af, 
the following relation can be obtained from equations (16) 
and (20): 

{/' " [>]*+,/*+.,} = (/ " RN)(M - SK MPC ){y sp - [y) k/k } 
+ (21) 

As the model used by the controller can differ from the real 
plant, this should be taken into account for obtaining the 
closed loop transfer function of the process. The real plant is 
assumed to follow the equations given below: 

[>]*+,«+. = M[y] klk + SAu k (22) 

where M and S are similar in form to M and S in equation 
(15). 

Note that for the real process there is no need to predict n r 
future instants. It is only necessary to evaluate h r future 
instants. Where h r is such that: 

h r T > max 0 :i 

So if the control actions are taken into account the relation 
given in equation (23) and the closed loop transfer function 
of equation (24) can be obtained. 

{y sp - } = [(/- RN)(M - SK MPC ) 

-RNSK MPC ]{y*-[y] k/k } 

+ RNM{y sp — [>']*/*} 

RNM 

M 



(23) 
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(24) 



where 

H = (I 



RN){M - SK, lf >c) - RNSK AtFC 

Disturbances can be easily incorporated in equation (24). 
This is made so as to evaluate the performance of the control 
strategy. One popular approach often used in the literature is 
to consider an input disturbance in the plant as was done by 
Morari and Lee . So, finally: 



H RNM 
~SK MPC M 







- M*/ 






-9k 



[du k ] 



(25) 



where du k is the input incremental disturbance. 

The transition matrix defined by equation (24) or (25) 
describes the closed loop poles and consequently the 
stability of the system. Also the eigenvectors of this 
matrix define the directions where the system is more 
sensitive to disturbances or where the error in the predicted 
output has a slower decay. Therefore, it is critical that the 



dimension of this matrix be reduced to the smallest possible 
value so that it is possible to perform an intensive 
verification of the system behaviour under different 
conditions. These should be taken not only in terms of 
changes in the parameters of the plant, but also in terms of 
the tuning of the controller. The location of the poles of the 
transition matrix defines clearly the behaviour of the system. 
This practically eliminates the need of extensive simulation 
of the process, as will be shown in section IV. 



Reducing the Order of the ROSSMPC — Final 
Formulation 

If it is wished to reduce the order of the system and 
preserve the stabilizing property of the long range MP 
controllers, the prediction model must be slightly modified. 
In the usual approach the predicted values of the system 
output are calculated at equally spaced sampling instants 
along the optimization horizon. With the model representa- 
tion of ROSSMPC the optimization horizon can be extended 
by including predictions of the system output at arbitrary 
future time instants. Or in other words a larger optimization 
horizon can be considered with fewer prediction points. As a 
consequence the state vector of the predictor becomes 
smaller. 

In the conventional approach the predicted output is given 
by: 

y k+J where j = 1,2 n h 

In the extended approach it is considered that the prediction 
is carried in the following instants j: 

;= i,2,. 

where 

n,T > max B tJ 



• >n n n t + n [y n t -f n 2 -f n k 



and n lt n 2 



>i k correspond to not necessarily equally 



spaced time instants and can be conveniently selected so 
as to take advantage of the stabilizing properties of the long 
prediction horizon. The idea of the method is to select a 
small horizon n t , complemented by another small set of 
sampling points n^n 2> n k . The complete optimization 
horizon rtjt+n, would be similar to the optimization horizon 
n r of the conventional MP controllers. The state variables 
are so now defined as: 

and the model coefficients of equation (15) are redefined as: 



M : 



*0 




0 


0 


0 . 


. 0 


0 


0 


0 


0 


/„ . ' 


0 


0 . 


. 0 


0 


0 


0 


0 


0 . 


■■ 7 «. 


0 . 


. 0 


0 


0 


0 


0 


0 


0 


0 . 


. 0 




Yn, 


0 


0 


0 


0 


0 . 


. 0 






0 


0 


0 . 


0 


0 . 


. 0 




***** 


0 


0 


0 . 


0 


0 . 


. 0 




0 


.0 


0 


0 . 


0 


0 . 


. 0 


0 


D 



Trans IChemE, Vol 75, Part A, October 1997 



ROSSMPC: A NEW WAY OF REPRESENTING AND ANALYSING PREDICTIVE CONTROLLERS 



699 



R = 



s = 



(26) 



Co 
CDN p 

Observe that for the case of extrapolation \p nr ^„ /+n ,,..., 
yp n +nt are calculated using equation ( 1 0) at the future prediction 

instants Oh + 1)7", («, + «, + 1)T («, + n k + 1)7. 

Though it is not an usual approach to consider time 
instants that are not multiples of the sampling time, 
ROSSMPC can easily accommodate these, which can be 
particularly useful for interpolating the system output at 
time instants between the sampling instants. So this is 
another advantage of the ROSSMPC over the other existing 
MP representations. 



IV. APPLYING THE ROSSMPC TO SEVERAL WELL 
KNOWN EXAMPLES 

In the previous section the ROSSMPC representation was 
introduced. There it was shown how the existing linear MP 
controllers (DMC, MAC, OBMPC, etc) can be transformed 
into the ROSSMPC formulation. By doing so a significant 
reduction in the number of states can be obtained, allowing a 
better knowledge of what is really happening in the closed 
loop system. This is clearly shown in example L There it 
can be seen that using the ROSSMPC representation it can 
be understood why the performance of the DMC controller 
used by Lundstrom et al 1 was poorer than the one of the 
OBMPC controller. Also the ROSSMPC representation can 
be used to better design the DMC controller and so similar 
performances will be obtained using either one of these 
controllers, which is shown in example 2. In example 3 the 
Shell problem is considered and it is shown how the 



ROSSMPC can be used to reduce the number of states of the 
system. 

Example-1: ROSSMPC as an Analysis Tool 

This example is intended to show how ROSSMPC can be 
used to analyse the performance of MP controllers. Here the 
DMC and the OBMPC controllers are considered. 

The example is taken from Morari and Lee 17 . The process 
is given by the following transfer function and is assumed to 
be exactly known, i.e. no model uncertainty is considered in 
this example: 



G(s) = 



IQQg-' 
1005+ 1 



(27) 



In order to achieve a high closed loop bandwidth, Morari 
and Lee 17 use a sampling time of lmin. As the steady-state 
is reached after about 300 min there is the need of about 300 
coefficients in the step response model. That is, the 
approach of Li et at. 1 will need 300 states for representing 
the process model. Therefore, for the closed loop simulation 
altogether 600 states will be needed (300 for the process 
model and 300 for the real plant), which is obviously 
inconvenient for the purpose of design and analysis. The 
approach of Hovd et al. 9 would need at least 100 time steps 
for the truncated step response plus a first order system to 
approximate the tail of the step response. Thus, to represent 
the plant and model of the state predictor more than 200 
states would be needed which is still excessive to represent 
this simple system. The number of states obtained by these 
two approaches is the same for the DMC and OBMPC. 
Using the ROSSMPC approach proposed in this paper a 
drastic reduction in the size of the state space model is 
obtained. For instance, for the DMC controller with a 
prediction horizon n r = 4. the compact model is applied with 
only 1 1 states as follows: 4 states for model predictions, 2 
states for model parameters. 3 states to calculate the plant 
output and 2 states as the new artificial states. Therefore, a 
reduction of over 90% in the number of states is obtained 
when the ROSSMPC representation is used. Note also that 
the response obtained by the ROSSMPC model is the same 
as the one obtained when using the conventional formula- 
tions of DMC, OBMPC etc. as can be seen from the figures 
presented next. 

As a compact step response model can be obtained, 
ROSSMPC may be used as an analysis tool to explain the 
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Figure 1. Responses for a step disturbance in ihe input: ( ) DMC and ( ) OBMPC. 
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performance of any controller. Next the approach for the 
DMC and OBMPC controllers will be exemplified, as these 
were used by Lundstrom et at? It will also be considered 
that the process is subjected to input disturbance and the 
closed loop eigenvalues and eigenvectors of the system will 
be calculated for the ROSSMPC compact form. The DMC 
controller is simulated for the following parameters: 

7=1; m = 2; n r = 4; h r = 3; A = 0 

The disturbance is applied at the sampling time =5. Figure 1 
shows the closed-loop responses, which are approximately 
the same ones obtained by Lundstrom et al? From the figure, 
one might think that there is a clear superiority of OBMPC 
for this particular example. In fact in example 2 it will be seen 
that this may not be the case. 

If in this case the eigenvalues of the transition matrix of 
equation (24) are inspected it can be observed that DMC 
drives the system poles to zero as a dead beat controller. The 
exception is a single pole that remains at 0.990 which 
corresponds to a time constant of 100 min for a sampling 
period of 1 min. In other words, the feedback strategy 
adopted by DMC for disturbance correction does not move 
the predictor pole from its open loop position. The 
normalized eigenvector that corresponds to this pole is: 

U/m/c = [y* +1 -.- y k+A p p y k+i y k+2 y M p p] 

= [-.0141 - .0071 0 .0070 .0731 -.7102 

-.0142 -.0141 -.0139 0 - .0141] 7 " 

This vector defines the direction where the error calculated 
by equation (25) decays with a time constant equal to 
100 min. This is exactly what happens for the output of 
DMC in Figure 1. It can be observed that components 1 to 4 
and 9 to 1 1 of v DMCy corresponding respectively to model 
and plant states, are small compared to components 5 and 6 
which correspond to the states of model parameters. This 
means that the slow approach of the error vector to zero is 
performed with small errors in the model and plant states. 

One last comment is noteworthy. The apparent large error 
of DMC in Figure ( 1 ) is related to the extremely high gain of 
the system represented by equation (27). A better criterion 
of performance should be to compare the disturbance 
rejection of both approaches observing the open loop 
response of the system. For the input disturbance case, after 
30 min in open-loop the system output would be about 25, 



while with DMC the system output is 1.5. This means that 
the controller rejects 94% of the disturbance effect, which is 
a reasonable performance for any practical standard. On the 
other hand, OBMPC rejects more than 99.9% of the same 
disturbance after 30min. This superior performance can be 
explained if the closed loop poles of OBMPC are also 
calculated. In this case the poles of largest module 
corresponding to the OBMPC of Figure 1 are: 0.7462 and 
0.0882 ± 0.1946/. In this case the error in the critical 
direction decays with a time constant of 3.5 min, or the 
closed loop with OBMPC is about 30 times faster than DMC. 

Thus, by calculating the eigenvalues and eigenvectors of 
the system, the two controllers can be compared without the 
need to simulate the system. This is of great interest as the 
calculation of the eigenvalues is a relative easy task and so 
the procedure could be used to design a controller, which 
will be presented in example 2. 



Example 2— Using the ROSSMPC for Better Design of a 
DMC Controller 

Suppose now that based on the closed loop information 
presented in example 1, it is required to improve the design 
of the DMC for the system of example 1 given by equation 
(27). Since the problem of poor performance has been 
diagnosed as the location of the open loop pole of the model 
used by DMC, the alternative to improve the performance is 
to use an observer faster than the real plant. As an example, 
consider a DMC with the following model: 



where, G M is the model used by the controller. 

Note that the model of equation (28) does not mean that 
the process is not exactly known. It simply means that a 
different model is chosen on purpose for the controller so as 
this is able to faster predict errors due to a disturbance in the 
inputs. If the process was not exactly known, the model (28) 
might not be the best choice. 

In this example again as done by Lundstrom et al? the 
process is assumed to be exactly known. Uncertainties in the 
process will be considered in sections V. and VI. 

The model of equation (28) was chosen so that there 
would be a time constant of the observer equal to one tenth 
of the plant time constant. For this DMC the closed loop 




0 10 20 30 40 50 "0 io 20 ST 40 50 
*"»(™> r.me(m.n> 
Figure 2. DMC responses for input disturbance:. ( ) with model of equation (28) and ( ) with ideal model. 
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model: 



12 14 16 IB 20 

Optimization horizon (ni) 

Figure 3. Example 3. Closed-loop spectral radius of conventional 
DMC( ) and ROSSMPC ( > with the extended horizon approach. 

poles become 0.8818, 0.2836 and -0.1754. The conse- 
quence of this is that the dominant time constant of the 
closed loop is about 8 min instead of 100 min. The responses 
of the system for the model of equation (28) and for the ideal 
case (model of equation (27)) are shown in Figure (2). 

It is clear from examples 1 and 2 that the analysis of 
closed loop poles is feasible with the proposed state space 
model and becomes a useful tool in the design of MPC. The 
position of these poles in the complex plane is a substitute 
for extensive closed loop simulations. 

Example 3: ROSSMPC— A More Convenient 
Representation 

The problem considered here is a part of the Shell 
problem (Prett and Garcia 10 ) and is given by the following 



£1- 



60s + 1 
4A2e' 2U 



5.5Se~ 21s 
50s+ 1 
7.20 



445+1 195+1 



[:;] 



The large difference between the time delays of the system 
outputs (0 and 28) makes DMC unstable for small 
optimization horizons. This is so, as for small prediction 
horizons the error on variable y 2 will be more significant 
since there is no dead time for it (and variable n 2 ). To reduce 
this effect the optimization horizon must be increased and 
consequently the dimension of the system will be enlarged. 
With the ROSSMPC approach fast and slow dynamics can 
be balanced without increasing significantly the size of the 
system. 

Figure 3 shows that the DMC with 7=5, m= 1, A = 0 
and r =12, is unstable for values of n r up to 15. In this 
figure is also plotted the spectral radius of the closed-loop 
transition matrix for the same controller, using the 
ROSSMPC representation extended with a single point 
«, + «, =25 and weighting the error at this point with 

=I °- ^ „ . kl , 

Observe that for this case the controller is stable tor any n t 

starting from n { = 1 which is the minimum value necessary 

to accommodate a time delay of 28 min. 

Figure 4 compares the responses of the system with the 

conventional DMC with n r =25 and the DMC under the 

ROSSMPC formulation with n, = 7 and n t + n, = 25. It can 

be seen that the performance is almost the same. In the 

conventional DMC approach, the transition state matrix of 

the closed loop would be of dimension 120 x 120= 14400, 




so 100 

Time (min) 




Figure 4. Example 3. DMC ( ) and ROSSMPC ( ) response for a change on the set points of the system outputs. 
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while the ROSSMPC approach with the compact system 
representation needs a transition matrix of order 
24x24 = 576. The reduction on the size of the system is 
enormous (96%). 

In this example it was shown that with the ROSSMPC 
representation of a controller can be obtained without any 
loss of performance! 



V. ROBUST TUNING OF MPC BY CONSTRAINING 
THE CLOSED LOOP POLES 

In the previous section, it was shown that a poor 
performance of MP controllers may be explained through 
the location of some closed loop poles. This is e.g. the case 
of the DMC controller of example I . For this example the 
poles were inherited from the open loop poles of the system. 
It might be thought that by looking at the poles of the system 
a controller could be designed in an adequate manner. This 
might be true when the model of the process is perfectly 
known. This is, however, not the general case. So a natural 
question is what will happen if there are uncertainties in the 
parameters of the plant model? 

Besides reducing the order of the closed loop system, 
ROSSMPC possesses some useful features that help in the 
analysis of the controller. One of these features is that all the 
model parameters, including the time delay, appear 
explicitly in the state-space formulation of the system. 
The other advantage is that each group of parameters (gains, 
time constants and time delays) appears linearly in the 
closed loop representation of the system. This fact will be 
used in this section and a tool presented of robustly 
analysing the stability and performance of the MP 
controllers. 

Remember that if there are errors in the identified model, 
these will only appear in the prediction model in the 
matrices M and S. So it can be assumed that the uncertainty 
in the matrices can be bounded as: 

S = S(I H~ A s ) 

M = M(I + A M ) 

When these errors are used, the transition matrix of equation 
(24) can be rewritten as: 



(/ - RN)(M - SK MPC ) - RNSU + A 5 )K MPC RNM(I 
-S{! + A 5 )K„ PC 
= A + EA A F 



/VA/(/+ A M )1 
MU + A„) J 



(29) 



where, 



A = 



(/ - RN){M - SK stPC ) - RNSK MPC RNM 



' —RNS 


RNM' 




~A 5 






; a, - 




-S 


M 







F = 



(30) 



no model uncertainties. 

x(k+\) = Ax{k) x(0) = x 0 (31) 

The eigenvalues of matrix A are supposed to be inside a 
specified closed contour BC of the complex plane. For 
stability reasons this contour must be inside the unit circle. 
Also note that an output multiplicative error is considered 
here as an example of how to robustly design a MPC by 
means of the ROSSMPC formulation. Any" other model 
error could be used. As will be shown later, uncertainties in 
the parameters of the model used by ROSSMPC can be 
easily put in the form of equation (29). 

The advantage of having the disturbed system being 
described by equation (29) is the possibility of using the 
stability radius concept of Hinrichson and Pritchard 11 .^This 
serves as a guide in the analysis of the closed loop stability 
and performance of the MPC described by the ROSSMPC 
formulation. Their result is summarized next. 

In equation (29) A had a particular structure, namely it 
was block diagonal. Let now A A be a full matrix. By doing 
the real stability radius of the transition matrix A can be 
obtained for a specified contour dC as given in Hinrichsen 
and Pritchard 1 1 as: 

r = inf {o(A A )\A A e R" 1 ™ 1 and A + EA A F 

has an eigenvalue at 9C} 

where a(A A ) is the largest singular value of A A and where ?tj 
n 2 are related to the structure defined by E and F. This 
definition is equivalent to 

r = inf inf{b(A A )\A A G /T tn - 

and det(sl -A- EA A F) = 0} 

It is straightforward to show that a manipulation of the 
determinant leads to 

r = infinf{a(A A )\A ii ^R n ' xn ~ 

*€dC 

and det[l - A A F(sI - A)~ l E) = 0} (32) 

Equation (32) can be so used to calculate the stability radius. 
Its form is, however, not so simple. Here the results of Li et 
al?~ become important, as they show a better way of how to 
calculate the real stability radius of real matrices, which is 
the case of the ROSSMPC formulation. 

Let <t>=Re<f>+jIm$=F(sI -A)~ ] E. Then Li et al xz 
show that for this definition the real stability radius can be 
calculated as: 



r s = < inf a, 



']}' 



(33) 



Note that matrix A corresponds to the transition matrix of an 
autonomous system described by equation (3 1) if there were 



Re4> —ylm<t> ' 
Im<f>/y Re<f> 

where o 2 is the second largest singular value of the matrix 4>. 
r=inf{r s ) (34) 

where, r s can be considered a local stability radius 
associated with the point s = x , + xj on the contour dC. 

Before making other considerations, it can now be shown 
how the above equations (33) and (34) can be used for the 
robust analysis of a given MP controller. 

In the design stage of the MPC, the usual available 
information consists of the nominal model and uncertainty 
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bounds on its parameters obtained through extensive plant 
data. The uncertainty could be then limited by a maximum 
value, e.g. a(A A ). The design problem would then be the one 
of selecting tuning parameters such that the system would 
remain stable for all possible plant models. Here plant 
models refer to the fact that when there are uncertainties 
several models could apply to the process depending on the 
actual operating conditions. 

When the uncertainties in the plant model can be limited 
the design procedure is given as: 

select the tuning parameters so that r > a(A A ) (PI) 

There are, however, situations where it is not known how 
large the uncertainties can be. For these cases the design 
procedure must be slightly modified. Here tuning para- 
meters must be selected that will enforce stability for the 
largest possible set of plant models, which can be viewed as 
the following procedure: 

mat r (P2) 

tuning parameters 

Thus the stability radius can play a major role in the 
robust stability design of a controller. Nonetheless some 
conservatism can be introduced if the stability radius is not 
adequately calculated. This will e.g. be the case when no 
structural information is taken into account. 

This is what happens with equation (33). There no 
structural information about the model error is considered 
besides that given by matrices E and F. Consequently, if 
equation (33) is used to calculate the stability radius, a 
conservative result is obtained as A A is considered to be a 
full matrix. So it is desirable to calculate the stability radius 
in a less conservative manner. Note that under the 
ROSSMPC formulation it was shown that if an output 
multiplicative error is considered A A becomes diagonal. For 
this case it will be shown next how to modify equation (33). 

For any diagonal matrix Q, the following identity is valid: 

det(l - A A 4>) = detU - QA A <f>Q~ l ) 
= det{l - A A Q<pQ~ l ) 
Making use of this identity, equation (33) can be rewritten 



r s = < inf < inf a 2 

1 0 1>G(0.1| 



Re4> —ylm4> 
Jm4>/y Re4> 



(35) 



where, 4> = n<t>il' ] - 

For the ROSSMPC, equation (35) can be used instead of 
equation (33) which gives better results. It remains now to 
show how errors in the parameters of the models can be 
expressed in the form of equation (30) so that equation (35) 
might be used. This will be done in section VI for the cases 
where there is model mismatch in the static gains and in the 
dead-time of the process. 



VI. APPLYING THE ROBUST DESIGN PROCEDURE 
TO SOME EXAMPLES 

In section V, the way the stability radius can be calculated 
was represented for the special model error structure 




•05 0 05 1 

Real pan of point on the unit circle 

Figure 5. Example 4. Local stability radius for ROSSMPC. 

presented in equations (29) and (30). It was also claimed 
that the ROSSMPC formulation allows the representation of 
the real system in this form. It was also shown that the 
robust tuning procedure of a MP controller under the 
ROSSMPC formulation can be stated as procedure (PI) or 
(P2). Here the aforementioned topics are illustrated in 
examples 4 and 5. 

In example 4 the following topics are covered: 

• there will be uncertainties in the static gains of the model: 
it will be shown how to obtain the structure of equation (30) 

• some tuning parameters of the ROSSMPC controller will 
be selected and how the stability radius is used to analyse 
the robustness of the closed loop system shown. First the 
stability case is considered and later it is shown how the 
performance can be analysed. 

• once the analysis tool is described procedure (P2) is 
applied to show how the parameters of the ROSSMPC 
controller can be tuned to achieve robust performance and 
stability. 

In example 5 the following aspects are considered: 

• there will be uncertainties in the time delay of the system: 
it will be shown how to obtain the structure of equation (30) 

• it is illustrated how the QDMC controller used by 
Zafiriou 18 and Lee & Yu 19 can be better analysed 

• once the behavior of the QDMC controller is understood, 
the effect of the mismatch in the time delay in the stability 
and performance of the system is shown. 




Figure 6. Example 4. Closed loop poles (+) for a plant with ft, = 0.271 1. 
k 7 = -0.27 1 1 and with ideal model (o). 
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Example 4 — Robust Analysis of the Shell Problem — 
Uncertainties in the Static Gains of the Model 

The process of example 3 is again considered, except that 
here it is assumed that there are uncertainties in the static 
gains as proposed by Prett and Garcia 10 : 

1 .77( 1 + it, )<r 28j 5.58( 1 + k 2 )e~ 27s ' 



60s + 1 
4.42(1 + k { )e~ 22s 



50s + 1 
7.20(1 +k 2 ) 



44s + 1 195+ 1 

For this special case it is noted that A A , = 0 and A s is a 
diagonal matrix given by: 



A, = 



As 



As a consequence, in the ROSSMPC structure E and F can 
be further reduced to: 



-RNS 



F=[K»pc 0] 



In the analysis that follows ROSSMPC will be applied to 
this system with the following tuning parameters: 

m = 2, n, = 7, 



T = 5, 



= 25, 



10 and A = 0 



(36) 



and the maximum values of the gain uncertainties ki and k 2 
will be calculated where the system still remains stable. 

Suppose further that any point of the unit circle can be 
represented by x, + x t i. Figure 5 shows the plot of r s 
calculated with equation (35) as a function of x r . 

From Figure 5 it is clear that the optimization problem 
defined in equation (34) is not convex as the function that 
relates r s and x r shows several local minimum points. The 
global minimum has to be selected by inspection, but this is 
not a major problem since the search is unidirectional and 
bounded. Then, it is concluded that for the selected 
controller, the stability radius of the system is r= 0.2711, 
corresponding to the point on the unit circle defined by 
x r = 0.810. This means that with the adopted tuning 
parameters, the system under ROSSMPC control is stable 
as long as the uncertainties obey the following condition: 
d(A s ) < 0.2711. 

These results can be interpreted by examining the poles of 
the closed loop system. Figure 6 shows the dominant poles 




-0.5 0 0.5 

Real pan of point on ecnicur delta C 

Figure S. Example 4. Performance radius for contour 8C, . 



of the closed loop system. There two situations are 
considered: 



• no model uncertainties 

• model uncertainties defined 
As = -0.2711. 



by *, = 0.2711 and 



It is observed that for the latter case there are poles at 
0.810 ± 0.587; which are exactly on the unit circle. Any 
further increase in the model uncertainty will push these 
poles outside the unit circle and the system with the selected 
ROSSMPC parameters will become unstable. 

Up to here only stability considerations were discussed. 
Another very important feature of any controller is its 
performance and moreover its robust performance. Here, 
the way the stability radius can be used to enforce robust 
performance is discussed. 

Performance criteria can be established by means of the 
suitable location of the closed loop poles on the complex 
plane. The ideal situation would be the case where all the 
poles of the closed loop system were at the origin of the 
complex plane. In this case the system output would exactly 
follow the set-point. As the ideal situation is not usually 
feasible, there must be a performance criterion. Here it is 
considered that a process has good performance if none of 
the following situations occur: 

• underdamped oscillatory systems with a small damping 
ratio. This corresponds to the case where the response of the 
system will oscillate for a long time before it finally 
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Figure 7. Example 4. Closed loop poles for ideal model (o) and contour dC 
( ) that define a performance region on the complex plane. 
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Figure 9. Example 4. Relation between the performance radius and input 
weight A. 
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Figure JO. Example 4. System outputs for a set-point change with A = 0 
( ) and A = 0.7 ( ). 

stabilizes. This kind of response is associated with poles 
with large negative real part or large imaginary part. 
Consequently this kind of response can be prevented by 
further constraining the poles of the closed loop into a 
suitable sub-region of the unit circle. 

• sluggish or overdamped systems with a large dominant 
time constant. This kind of response is related to the 
presence of closed loop poles on the real axis close to (1,0). 

The robust performance design procedure could thus be 
viewed as the one of selecting a contour on the complex 
plane and of finding tuning parameters, that would produce 
a stability radius larger than the model uncertainty. 

Next the way to evaluate performance requirements for 
this example will be discussed. The contour represented in 
Figure 7. named 3C,, are selected as our performance 
criterion. It will limit a certain region of the complex plane, 
which will have the characteristics described next. 

On the left half side of the complex plane, contour dC x is 
a circle with radius r, = 0.65 and on the right half side it is 
an ellipsis with axes r, = 0.65 and r 2 = 0.95. The region 
limited by these contours enforces the following character- 
istics on the system: 

• the closed loop system will have a damping ratio of at 
least 2.5 for any plant model. Also the closed loop will never 
have a time constant larger than 78 min. 

For simplicity of notation, the stability radius related to 
this contour will be called the performance radius. Figure 8 
shows the local performance radius of ROSSMPC with the 
parameters given as in equation (36), as a function of x r for 
contour 3C,. By inspecting this curve it is observed that the 
performance radius has the value r, = 0.0818 at x r = 



-0.1885. This means that for any uncertainty A 5 for 
which a(A s )^0.0818, all the closed loop poles of the 
system will stay inside the region defined by contour BC X . 

In equation (36) A was considered to be equal zero. For 
this case it was seen that the system is to be stable. Except 
that the bound on a(A s ) will not be large. So a natural 
question would be how to better select A so that a larger 
limit on the uncertainty is obtained. The procedure that must 
be used to answer this question is the one given in (P2). That 
is what will be done next. 

For simplicity it is assumed that the two inputs have the 
same weight so that A will have the same diagonal entries. 
So whenever a value is assumed for A, it means that its 
diagonal elements have this value. 

Figure 9 shows the effect of input weight A on the 
performance radius. It is observed that the performance 
radius can be maximized by properly selecting the input 
weight. For this example the performance radius has a 
maximum value at A = 0.7, where the performance radius is 
close to rj = 0.17. 

This result means the following: 'the maximum error in 
the static gains must obey the following condition 
a(A 5 )=0.17 if the system is to attain the performance 
criterion Note that the controller with A = 0, does not 
attend the adopted performance criteria as a(A s ) = 
0.17 > 0.0818. This does not mean that it will be unstable, 
only that the performance will not be acceptable for larger 
errors in the static gains. To better illustrate this point 
consider next that the uncertainties in the static gains are 
defined by *, = * 2 = 0.17. Figure 10 shows the response of 
the system with model mismatch for a set-point change, the 
system being controlled by the ROSSMPC controller tuned 
with A = 0 and A = 0.7. The responses of output >•, are quite 
similar for both input weights, but for output y 2 the system 
response does not have the required damping ratio when 
A = 0. 

Example 5— Effect of Uncertainties in the Time Delay 

Consider the SISO system used by Zafirou 18 and Lee and 
Yu 19 to study the loss of robustness of MPC in the presence 
of uncertainty in the time delay, when there is hard 
constraint in the system output. 
1 

G(s)- 



-0.1 5j 



5+1 



In this case the step response of the nominal system is: 

J( r)=l-i.^- 015) 
While the true plant model is represented by: 

where A0 is the time delay mismatch between the true plant 
model and the model used by ROSSMPC. With this 
assumption the plant step response can be written as 
follows: 

5 = 5+ VAS 
where 

V= [_,-e™» 
AO = e~" - 1 
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Figure 1 /. Example 5. Stability radius ( ) and performance radius 

for contour dC, ( ), 

Applying these definitions in equation (29) we can obtain 
the following relations: 



-RNV RNW X 
-V W, 





" A6 0 " 




K\t/>c 0 




0 AO 




0 W 2 _ 



(37) 



where 
and 



0 0 ... 0 e -[(n t + l)T-0.\5] q q 



\V 2 = [0 ... 0 1] 

It is easy to see that matrix E in equation (37) is rank 
deficient and consequently matrix </> = F(sl - A)~ l E is 
singular for any point s on the chosen contour dC. As 
det(l — A A <f>) = 0 it can be shown that the stability radius 
associated with point s can be calculated by: 

r s = |AS| = I !/(*,., +^)| (38) 

where <f> u and <f> 22 are the diagonal elements of matrix <f>. 

Zafiriou and Lee and Yu 19 used the QDMC as the 
controller for this system. They have also shown by 
simulation that when AS = 0.15, the closed loop is unstable 
if a hard constraint in the system output becomes active. 
Here this will be analysed by means of ROSSMPC. The 
nature of critical values of the delay mismatch for both 




Figure 12. Example 5. Closed loop poles for &$ — 0 (o); A$ = -0.0513 0c) 
and A0 = -0.0260 (+). 
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Figure 13. Example 5. System input and output for a set-point change* 
AO = 0 ( ); &d = -0.0513 ( ) and = -0.0260 ( -). 

robust stability and robust performance will be studied. 

Robust performance will be defined in the same way as in 
the previous example, i.e., the system will have a good 
performance if the closed loop poles are within the contour 
BCi defined in example 4. 

The authors cited used the QDMC controller with hard 
constraint in the output. So first it must be ascertained how 
to formulate this controller using the ROSSMPC. The 
situation corresponding to having in the QDMC the 
constraints in the output active, is the one of the dead-beat 
tuning of ROSSMPC. With the least squares control 
objective function, this is obtained when the input 
suppression factor equals zero and the control horizon 
equals the number of prediction instants beyond the time 
delay in the optimization horizon. If there are no constraints 
in the input, it is possible to obtain control actions, that will 
eliminate the error between the predicted trajectory of the 
system output and its set-point. 

So, consider now the ROSSMPC with the following 
tuning parameters: 

7* = 0.1, A = 0, m ~ 5 and n t = m + 2 

Results are independent of the control horizon m, as long as 
the optimization horizon is equal to the control horizon plus 
the time delay. The stability radius of the system calculated 
by equation (38) is shown in Figure 1 1 as a function of the 
real part of the point s on the unit circle. The minimum of r s 
is equal to 0.0488 and is obtained for x, = -1. This means 
that the controller can be unstable if 

|A9| = \e** - 1 1 > 0,0488 

The critical case corresponds to AO = -0.0513 where 



Trans IChemE, Vol 75, Part A, October 1997 



ROSSMPC: A NEW WAY OF REPRESENTING AND ANALYSING PREDICTIVE CONTROLLERS 



707 



one of the closed loop poles lies at (—1,0), as shown in 
Figure 12. This explains why QDMC is unstable when 
AS = 0.15, as shown by Zafiriou 18 . The value of the 
uncertainty used by that author in his simulations was much 
larger than the allowable value that can be obtained by 
inspecting the stability radius of the system. 

Analogously to the procedure followed in example 4, a 
performance criterion can also be specified for this system, 
by constraining the poles of the closed loop system into a 
suitable closed contour of the complex plane. Adopting the 
performance contour dC t , the critical mismatch in the time 
delay becomes AS = —0.0260. The dominant pole of the 
closed loop system lies at (—0.65,0). This is also shown in 
Figures 1 1 and 12. 

To better clarify these results, Figure 13 shows the 
response of the system for a set-point change for the 
following cases: 

• no model mismatch 

• a model mismatch given by A6 = —0.0513 

• a model mismatch given by AS = —0.0260 

This example shows that ROSSMPC has a significant 
advantage over the other MPC representations when robust 
stability and performance of the system is the issue for time 
delay uncertainty. The model representation of ROSSMPC 
allows a simple and straight approach to solving the 
problem. 

VII. CONCLUSIONS 

In this paper a state space representation has been 
presented for the class of model predictive controllers, 
named ROSSMPC. This representation can be used for any 
existing MP controller. For example, ROSSMPC was 
applied for the DMC, QDMC and OBMPC controllers. So 
ROSSMPC can be also regarded as a way of describing 
existing MP controllers. This is of great interest as several 
industrial applications make use of the DMC controller. 
While preserving the properties of the existing controllers 
on the one side, ROSSMPC permits us to better explore and 
comprehend the controller chosen. This is so as the 
ROSSMPC is based on a step-response model and preserves 
the parameters of the model in an explicit manner in the 
state space representation. Moreover the number of states is 
reduced. As a consequence the eigenvalues or poles of the 
closed-loop transfer function can be evaluated and a better 
comprehension of the behaviour of the closed-loop system 
is acquired without the need of extensive simulation. This 
also makes the procedure of tuning the controller simpler. 
The parameters of the controller can be so chosen as to make 
the poles of the transfer function lie in desired regions of 
the complex plane. For the ROSSMPC representation, the 
prediction of the behavior of the system in the future can be 
arbitrarily made and there is also no need of using equally 
spaced time instants. For processes where there is the 
necessity to use large control horizons, this is of great 
advantage. Stability can be ensured and at the same time the 
number of states remains small. ROSSMPC can be used 
both as an analysis tool and as a way of modelling and 
designing a controller. For the first case a given MP 
controller can be better designed. 

For example, in this paper it was shown how the model 
used by the DMC controller can be modified in order to 



enhance its performance. That is, it was shown that 
sometimes introducing known errors in the model of the 
process can make the controller respond faster to dis- 
turbances. This is an interesting result and shows that the 
model itself can be regarded as another parameter that can 
be selected to enhance the performance of control systems. 
As the ROSSMPC model possesses fewer parameters than 
other MP models, it can be used as a controller and so the 
computer work of calculating control actions is reduced. 
This is important for the cases where there is interest in 
having one or very few controllers for the whole unit. 
ROSSMPC has a structure that is also suitable for the use of 
robust analysis tools. 

In this paper it was shown how the stability radius can be 
used for allowing the controller to work well and maintain 
robust stability and performance. This was, for example, 
accomplished by limiting the closed loop poles inside 
suitable regions of the complex plane. 



NOMENCLATURE 



A 


transition matrix for the closed loop with nominal model 


c 


parameter of the continuous step response function 


C 


matrix of parameters of the step response function 


d 


coefficient for updating the prediction function parameters. 




equation (12) 


D 


matrix of coefficients for updating the prediction function 




parameters 


E,F 


matrices defining the structure of the uncertainty, equation (30) 


G 


process transfer function 


f 


function in the continuous step response, equation (3) 




MPC feedback gain 


m 


control horizon 


M 


state matrix in the prediction model, equation (4) 


M 


state matrix of the plant real model 


n h 


step or impulse response model horizon 


*r 


optimization horizon 


"r 


first part of the prediction horizon of ROSSMPC 


h t 


prediction horizon of the plant in the reduced order model 


n u 


number of inputs 


n y 


number of outputs 


N 


matrix in the error correction equation (5) 




parameters of the output prediction function, equation (S) 


R 


disturbance correction coefficient for DMC, equation (5) 


S 


step response coefficient of the model in MPC 


S 


step response coefficient of the plant 


S 


step response coefficient matrix of the prediction model, 




equation (4) 


s 


plant step response coefficient matrix, equation (22) 


s 


dynamic matrix of MPC 


t 


time 


T 


sampling period 


u 


manipulated input 


u 


Eigenvector of the closed loop transition matrix 


y 


system output 


Creek symbols 


r 


output weighting matrix, equation (19) 


A 


input weighting matrix, equation (19) 


e 


time delay 




matrix in the output prediction, equation (10) 




singular value 




structured singular value 


Abbreviations 


DMC 


dynamic matrix controller 


MAC 


model algorithmic control 


MP 


model predictive 


MPC 


model predictive controller 


MPHC 


model heuristic predictive control 


OBMPC 


observer based model predictive control 


QDMC 


quadratic dynamic matrix control 


ROSSMPC 


reduced order state space model predictive controller 
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